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Abstract — Based on autoregressive models of neural network 
identification the equation of one-step and multi-step predictor 
variables showed on the basis of NNARX, NNOE models. 
Research of the influence of the depth of regression on the 
quality of the identification and on the speed of convergence was 
performed. The results of identification modelling on real 
performance of beneficiation process are given below. 

Index Terms — identification; NNARX; NNOE; beneficiation 
process, magnetite quartzite. 


I. INTRODUCTION 

Opening question regarding the identification of 
technological process (TP) parameters of magnetite quartzite 
(iron ore) beneficiation process in industrial environments for 
next intelligent regulator synthesis for all production stages 
coordinated control. The relevance and the statement of this 
problem are given in the articles of the authors [ 1 -2]. 

Previous researches results showed that neural networks 
approach can be chosen for identification of technological 
parameters of beneficiation process under conditions of ore 
concentration plants. Required quality and convergence of the 
model can be reached, in case of correct selection of 
modelling structure basis. There is also an open issue of some 
parameters of neural network models selection, which directly 
influence on the quality and speed of the system training 
(parametrization). Resulting from the described above, the 
purpose of current article is the definition of optimal time 
delay values of input and output signals for autoregressive 
models of NNARX, NNOE [2-3] type. 

n. Models of Identification and Problem Statement 

In accordance with [4-6] TP of beneficiation process 
considered as multidimensional discrete system of the form: 

Z p = {[w(£), v(&)], k =1, iv},(l) 

where P is order of system; u(k) is vector of input control 
signals; y(k) is vector of output signals of system; N is 
maximal quantity of discrete supervision readouts; k is 
discrete time. 

Autoregressive NNARX and NNOE models are mostly 
widespread now. According to [3], they are recommended for 


use in case of a limited set of input statistics (e.g. creating an 
automated process control system). 

Generally, the mathematical models (in other way, 
approximators, predictors or regressors) based on artificial 
neural networks are represented by nonlinear equation of the 
form: 


y(k + l) = NN 


y(k),y(k -1),..., y(k -/j), N 
u(k),u(k -l),...,u(k -l 2 ) 


( 2 ) 


where y is output signals vector of identificational model; 
NN(-) is resulting transformation of "input-output" type, 
which is performed by the neural network; 11 is regression 
depth (the number of delay signals) of feedback on the output 
of the model; 12 is input delay depth. 

This is one-step predictor or NNARX model (Neural 
Network based AutoRegressive exogenous signal). Such 
predictor allows forecast one step forward and can be 
successfully used for identification. 

To solve this problem, neural covariates models, which 
allow to perform multi-step forecasts (short-term predictors) 
[6-7] are quite successfully used too. 

NNOE (Neural Network Output Error) is one of such 
models, which defined by recursive equation of the form: 


y(k + [) = NN 


y(k), y(k - 1),..., y(k - /,), u(k), 


■ (3) 


yi(k -1),..., u(k —1 2 ) 

As model structure of approximator, architecture based on 
the recurrent network with memory is used (Figure 1). 
According to the assumptions, stated in [8-11], these 
particular neural network structures are quite effective for 
modelling of beneficiation TP. 

b P“ - Output 



Figure 1. Analyzed structures of neural networks 
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III. Technique of Researches 

For parameterization (training) of models types (2)-(3) 
certain amount of statistical information is needed, namely 
data templates of "input-output" type. These data templates 
represented as multitude: 
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i = l,M, (4) 

where z is input template; d is output template; M is 
specified number of templates. 

Hence, according to (1-4), vector of input signals of the 
system is defined as: 

r y(k), y{k - 1),..., y(k - 
y ii(k -l),...,u(k -l 2 ) 

and the resulting output vector d = y(k +1). 

Best regressor model must provide an acceptable training 
rate, generalizing properties, and on the other hand requires 
proper prediction accuracy. 

To assess the adequacy of the resulting neural network 
model it is reasonable to use expression: 

||y-;y||<£, 

where • is determined norm; e is measure of approximation 



inaccuracy. 

Traditional criterion here is average quadratic error of 
prediction (quadratic norm) of the form: 

i M i M 

MEE =- >'W ) 2 =777 .• 


t =1 


2 


To assess the depth of regression, method of Lipschitz 
coefficients is used, which described in [12]. Evaluation of 
delayed input and output signals quantity is made on the basis 
of such criteria: 


q (l) = 


W\k) 


\k =1 


i Ip 


, (5) 


where q {l) , q U) (k) are respectively current average 
Lipschitz coefficient value under delay depth in / -signals 
(considered that /j = l 2 = l ); p is the quantity of largest 
selected coefficients out of all possible "input-output" type 
combinations of pairs ( p = 0,0 IN 4- 0,02N ). 

Current Lipschitz coefficients values are calculated 
according to the formula: 


r ( 0 = y(kj)-y(kj) 

uik^-uikj) 


i * j. 


where • denotes the Euclidean norm, and the values of the 

regression input and output vector are determined based on 
equations (2) and (3). 

We perform simulations of TP identification for the magnetite 
quartzite beneficiation, using covariates of NNARX (2) and 
NNOE (3) types and perform selection of delay signals 
optimal quantity by the criterion (5). 

According to [1] during the progress of magnetite quartzite 
beneficiation process, under conditions of modern ore 
concentration plants it is necessary to control next set of 
factors: 

^p _ ^ ^0’ 2o’ Q-> d, P m , p K , p c , 1 

1 0 m ’0„’0K>r,r K >£,£ K ]’ 

where (X = \(X j }, i — \...N r is useful component content 


value in raw ore (in magnetite quartzite batch); N r is quantity 


of industrial species in the batch; ^ } is density of the 

respective species; g = \g : } is a proportion (ratio) of 

species in the batch; do is averaged ore coarseness before 
beneficiation process; Qo is ore feed at the first stage of 
beneficiation process; Q = )Q is 

processing at each stage of beneficiation process; N s is 
quantity of beneficiation process stages; d = \d - } is 
averaged ore coarseness of a product after each stage of 
beneficiation process; P m = \P mj } is content of solid 

material during degradation process; P K = \p /: ^ | is a density 
of pulp before classification process; p c = \P C - j is a 
density of pulp before magnetic separation process; 
Ppp = ^pp,) is content of useful component in industrial 


product after each stage; J3 x = \P XJ } i 


is losses of useful 


component in tailings; j3 k is quality of concentrate; 


r = 


W is 


is an output of useful component in final middling 

product; y K is an output of useful component in concentrate; 
£ = \p . | is extraction of useful component in final middling 

product; S K is extraction of useful component in concentrate. 


IV. Modelling Results and Conclusions 

Training of neural networks requires a statistical 
representative sample of TP. According to these 
requirements, during 2006-2016, authors collected statistics 
of basic technological performance indicators of all mining 
complexes concentration plants in the city of Kryvyi Rih, 
Ukraine (Tsentralnyi, Inguletskyi, Pivnichnyi, Pivdennyi, 
ArcelorMittal Kryvyi Rih mining complexes). The list of 
indicators was defined according to (6). Larther on, according 
to the method of passive experiment, preliminary data 
processing was done. The individual "emissions" (incorrect 
data) was removed. Then, according to the requirements of 
(4), the set of indicators for each plant was divided into two 
subsets: 

1) subset for training (parameterization) models (up to 2000 
templates); 

2) subset of the tests (verification of adequacy) (up to 500 
templates). 

Subsequent modelling was based on the described 
methodology by using Neuro Solution package. The recurrent 
network model was taken as basis with the usage of direct 
signals and elements with memory distribution technology 
(figure 1). Lor parameterization of models, 
Levenberg-Marquardt training method (LM - method) was 
used. 

Table 1 and figure 2 show main results of simulation the 
identification of magnetite quartzite TP beneficiation process 
and calculation of Lipschitz coefficients for different models. 
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Figure 2. Dependence of the Lipschitz coefficients values(g) on the 
depth of regression (/) for NNARX, NNOE types of models 


TABLE I. Quality and convergence indicators of 

MATHEMATICAL MODELS 


JV» 

Indicators and types of models 

Indicators name 

NNARX 

NNOE 

1 

MSE (maximal) 

2,62 

2,58 

2 

Training speed (in epochs) 

588 

893 

3 

Training method 

LM 

LM 
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V. CONCLUSIONS 

On the basis of modelling results, conclusions were drawn. 

1. For identification of TP of iron ore magnetic beneficiation 
process, usage of NNARX and NNOE type of neural network 
predictor models is sufficiently effective. 

2. In the process of identification with usage of NNARX type 
of model sufficient depth of delay is 2-3 signals (based on 
point of fracture position). 

3. Usage of NNOE type of model requires a marginally 
greater depth of regression (3-4 signals). 
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